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PREFACE 

This  is  an  interim  report  on  an  active  research  field.   It 
contains  a  comprehensive  coverage  of  those  factors  which  influence 
LSA  ignition,  provides  a  useful  introduction  to  the  mathematical 
treatment  of  LSA  wave  propagation  and  reaches  certain  conclusions 
about  LSA  wave  ignition. 

During  the  period  between  completion  of  the  manuscript  and 
publication  of  this  report,  investigation  of  these  problems  has 
continued  and  thus  many  of  the  speculations  reported  here  have 
since  been  confirmed. 


UNCLASSIFIED 


UNCLASSIFIED 

CONTENTS 

Page 

1.0  INTRODUCTION  AND  SYMBOLS  1 

1.1  SYMBOLS  3 

2.0  PLASMA  IGNITION  5 

2.1  IGNITION  MECHANISMS  6 

2.11  Atom  Emission  6 

2.12  Ion  Emission  8 

2.13  Electron  Emission  10 

2.14  Photon  Emission  12 

2.15  Thermomechanical  Response  12 

2.2  IGNITION  EXPERIMENTS  13 

2.21  Laser  Supported  Combustion  Wave  Ignition  14 

2.22  Laser  Supported  Detonation  Wave  Ignition  15 

Table  I  17 

Table  II  18 

Table  III  19 

2.23  Plasma  Emission  Experiments  20 

Figure  1  20A 

Figure  2  21A 

Figure  3  21B 

3.0   QUALITATIVE  SUMMARY  OF  THE  TRANSITION  PROCESS  23 

Figure  4  23A 

4.0  PROPAGATION  OF  LSA  WAVES  27 

4.1  HYDRODYNAMIC  EQUATIONS  28 
Figure  5  28A 

4.2  ONE  DIMENSIONAL  LSD  WAVES  33 


UNCLASSIFIED 


UNCLASSIFIED 

4.3  ONE  DIMENSIONAL  LSC  WAVES  36 

Figure  6A  36A 

Figure  6B  36B 

Figure  7  38A 

Figure  8  38B 
Figure  9                               .         .  39A 

Figure  10  40A 

4.4  TWO  DIMENSIONAL  LSA  WAVES  41 


5.0   SUMMARY  42 


UNCLASSIFIED 


UNCLASSIFIED 


LASER  ABSORPTION  WAVE  PHENOMENA 

Don  E.  Harrison,  Jr. 
Physics  and  Chemistry  Department 


1.0   INTRODUCTION  AND  SYMBOLS 

(U)  The  topics  discussed  here  appear  in  the  literature  under  a 
variety  of  headings.  Chief  attention  will  be  focused  upon  Laser  Supported 
Detonation  (LSD)  waves  and  Laser  Supported  Combustion  (LSC)  waves.   It  is 
convenient  to  group  these  two  phenomena  under  the  more  general  heading  of 
Laser  Supported  Absorption  (LSA)  waves,  which  includes  both  surface  init- 
iated and  gas  initiated  LSC  and  LSD  waves  as  well  as  plasmons  or  plasma- 
trons,  which  describe  a  related  phenomenon,  a  plasma  which  does  not 
detach  from  the  target  surface.  The  best  available  summary  of  current 
thinking  on  these  topics  is  Volume  V  of  the  Proceedings  of  the  October 
1973  LASER  EFFECTS/HARDENING  CONFERENCE  in  Monterey,  California1. 
Warning :  Understanding  has  increased,  and  many  individual  positions  on 
controversial  topics  have  been  modified  since  its  publication! 

(U)  The  introductory  article  in  reference  1  by  P.  E.  Nielsen2  is  a 
valuable  overview  of  many  of  the  topics  discussed  here.  The  entire  field 
of  study  depends  to  a  large  extent  upon  the  initial  analysis  by  Raizer3, 
published  in  1965.  A  useful  discussion  of  the  relationship  between  the 
various  phenomena  is  contained  in  a  report  by  Nielsen  and  Canavan*4  which 
is  unpublished. 

(U)  This  article  considers  a  family  of  effects  which  appear  when 
laser  radiation  with  irradiance  of  the  order  of  105  to  108  watts/cm2  is 
incident  upon  a  solid  target  in  air  at  pressures  above  approximately 
0.1  atm.  Under  these  conditions  a  plasma  forms  at  the  target  surface 
and  then  may  move  back  along  the  laser  beam  if  the  irradiance  is  suffic- 
iently large.   If  the  plasma  expansion  along  the  beam  is  subsonic,  the 
phenomenon  is  called  an  LSC  wave;  if  it  is  supersonic,  the  wave  is  an 
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LSD  wave.  The  general  propagation  characteristics  of  the  two  types  of 
waves  have  been  known  for  some  time,  but  significant  aspects  of  their 
propagation  and  initiation  are  still  under  active  investigation.  For 
present  purposes,  initiation  has  been  divided  into  two  steps:  ignition 
in  which  the  first  few  ions  and/or  electrons  are  produced,  and  trans- 
ition which  carries  the  initial  ionization  into  a  full  blown  plasma  which 
is  propagating  as  an  LSA  wave. 

Qj)  LSA  ignition  occurs  at  irradiance  levels  which  are  two  orders 
of  magnitude  lower  than  those  normally  associated  with  gas  breakdown 
and  spark  formation,  the  gas  initiated  LSA  phenomena.  The  laser-target 
interaction  is  poorly  understood;  so  a  large  number  of  ignition  mechanisms 
have  been  proposed.  Hopefully  this  review  can  decrease  the  size  of  this 
list  of  candidates. 

(U)  Once  ignited,  the  transition  to  full  LSA  wave  status  occurs  by 
well  understood  physical  processes.  Unfortunately  the  details  of  the 
transition  depend  upon  the  characteristics  of  the  laser  (particularly 
its  wavelength)  and  upon  the  interaction  of  a  number  of  competing  pro- 
cesses. The  role  of  each  process  can  be  identified,  but  its  relative 
contribution  to  the  transition  phenomenon  is  uncertain.  As  a  result,  no 
detailed  theory  of  transition  for  a  particular  target-laser-atmosphere 
system  exists,  even  though  transition  times  can  be  estimated  with  some 
confidence. 

(U)  This  article  surveys  the  propagation  of  LSA  waves  so  that  the 
differences  between  LSC  and  LSD  waves  can  be  clearly  established.   It 
also  provides  background  information  about  hydrodynamics. 
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1 . 1  SYMBOLS 


The  symbols  listed  here  are  in  the  order  of  their  appearance  in 
the  article. 


X 

wavelength 

p 

density 

sub: 

scripts 

m 

metal 

V 

vapor 

1 

ambient  air 

2 

shocked  air 

v  speed  of  vaporizing  surface 

G  irradiance 

L  heat  of  vaporization  of  metal 

C  specific  heat  of  metal 
m      r 

T  temperature 

j  ion  current 

j  electron  current 

0+  positive  ion  surface  work  function 

4>  electron  work  function 

I  ionization  energy  of  an  atom 

Z  atomic  sublimation  energy  from  metal  surface 

j  vaporization  current  of  neutral  atoms  from  a  target  surface 

k  Boltzmann's  constant 

v  frequency  of  radiation 

T  lifetime  of  an  electron  heating  cycle 

e  electron  charge 

N  electron  density 
e  ' 

3  UNCLASSIFIED 


UNCLASSIFIED 


N  critical  electron  density  for  opaque  plasma 

P  pressure 

V  specific  volume,  that  is,  volume  per  unit  mass  V  =  1/p 

e  specific  internal  energy 

u  fluid  flow  velocity  in  LAB  system 

u*  fluid  flow  velocity  in  COM  system 

h  specific  enthalpy 

t  volume  element 

D  shock  front  propagation  velocity 

* 

q  =  G/p2u2 


j 

*         * 

=  u2p2  =  uipi 

c 
p 

Oh/3T)p 

c 

V 

(3e/3T)v 

a 

C  /C 
p   V 

c 

speed  of  sound 

K 

absorption  coefficient 

0 

thermal  conductivity 
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2.0  PLASMA  IGNITION 

(U)  At  least  ten  ignition  mechanisms  have  been  proposed  since  1971. 
In  some  cases  the  proposal  consists  of  little  more  than  a  name,  while  in 
others  a  relatively  complete  theory  has  been  developed.  Many  of  the 
ideas  were  developed  through  private  communications;  so  it  is  exceedingly 
difficult  to  correctly  attribute  them  to  their  originators.  The  situa- 
tion is  further  complicated  by  the  intermixing  of  ignition  and  transition 
processes  in  the  presentation  of  a  particular  mechanism.   Both  Nielsen2 
and  Walters6  attempt  to  give  proper  credit;  so  these  references  should 
be  consulted  by  the  historian. 

(U)  This  review  will  take  another  direction.  The  ignition  mechanisms 
will  be  grouped  by  type,  and  then  the  common  and  distinctive  features  can 
be  more  readily  identified.  The  five  types  are:  atom  emission,  ion  emis- 
sion, electron  emission,  photon  mechanisms  and  shock  wave  mechanisms. 
This  organization  will  facilitate  the  description  of  the  ignition  pro- 
cess. 

(U)  The  experimental  fact  is  that  a  plasma  forms  near  (within  500 
urn)  the  target  surface  within  a  very  short  time.  This  time  is  of  the 
order  of  50  nsec  for  LSD  waves  and  somewhat  larger  for  LSC  waves.  The 
plasma  contains  ions  from  both  the  target  and  from  the  atmosphere .   In 
all  probability  it  is  formed  by  ions  attributable  to  both  sources,  but 
there  is  some  evidence  that  atmospheric  ions  appear  first  in  the  LSD 
environment  while  target  ions  appear  first  in  the  LSC  environment.   In 
either  event,  the  transition  plasma  feeds  upon  both  sources,  and  the 
plasma  which  can  be  examined  spectroscopically  contains  ions  from  both. 
Because  of  the  nature  of  this  review  both  possibilities  will  be  discussed 
in  some  detail;  so  that  the  information  will  be  available  regardless  of 
the  ultimate  experimental  resolution  of  the  present  dilemma.  The  avail- 
able experimental  conclusions  can  then  be  discussed  in  context. 

(U)  The  ignition  .of  LSA  waves  is  closely  related  to  the  more  general 
problem  of  laser-induced  particle  emission  from  surfaces.  Many  experi- 
ments have  been  performed  with  ruby  lasers  (X  =  0.6943  ym  ■*  1.78  eV 
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photons)  which  are  quite  comparable  to  the  Nd-YAG/glass  (X  =  1.06  ym 
♦1.17  eV)  and  C02  (A  =  10.6  ym  ♦  0.117  eV)  experiments  reported  by  the 
LSA  specialists.  These  investigations  have  been  surveyed  by  Ready7, 
whose  book  should  be  consulted  for  further  details. 

(U)  As  a  general  principle,  ignition  effects  rooted  in  thermal 
processes  should  vary  little  between  ruby  and  CO2  wavelengths.  The 
thermionic  emission  of  electrons  must  fall  into  this  group.   In  contrast, 
photon  sensitive  processes  such  as  the  photoelectric  effect  would  hardly 
be  comparable  because  of  the  difference  in  photon  energy.   In  the  suc- 
ceeding sections  ruby  laser  results  will  be  included  when  appropriate. 

2.1   IGNITION  MECHANISMS 

2.11  (IT)  Atom  emission  from  the  target  is  one  plausible  source  for  the 
plasma  if  ionization  mechanisms  are  also  available.  The  atoms  may  be  of 
the  target  species,  from  a  surface  oxide  layer  (in  which  case  the  emitted 
species  may  be  a  molecule,  A10,  for  example),  from  surface  included 
material  (laid  down  during  polishing  and  finishing) ,  or  from  simple 
•dirt'  which  can  find  its  way  to  the  surface  in  untold  number  of  ways. 

(U)  There  will  be  a  characteristic  vaporization  temperature  for 
each  of  the  possible  species,  a  specific  heat  for  each,  and  a  latent 
heat  of  vaporization  for  each.  Momentum  conservation  guarantees  that 
the  emission  process  would  be  thermal,  which  implies  that  the  vaporiza- 
tion temperature  determines  the  velocity  distribution  of  the  emitted 
atoms.  The  mean  emission  speed  of  the  vapor  is  determined  by  the  con- 
servation of  mass  and  energy  and  the  parameters  of  the  source  and 
emitted  material .  . 

(U)  As  an  example,  consider  emission  from  a  perfectly  clean  metal- 
lic target.  Then  conservation  of  mass  requires  that: 

p  v  =  p  v  ; 
Km  m   *v  v  ' 

here  v  is  the  speed  with  which  the  vaporization  surface  penetrates  the 
metal,  and  conservation  of  energy  requires  that: 
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v  =G/(p  [L  +  CT]). 
m     v  mL  m   m  mJ 

The  vaporization  speed  is  difficult  to  determine,  so  use  a  relation 
proposed  by  Anisimov,  et.al.8: 

v  =  c  exp-(L  /CT). 
m    m   r  v  m  m  nr 

After  some  algebra  the  vapor  speed,  v  becomes 

v  =  AL   ,  (1) 

v     m  '  v  J 


with 


A  =  [2MP   L  c  /G)]*  . 


m  m  m 


The  emission  speed  is  of  critical  importance  because  it  scales  the  time 
available  before  the  atom  leaves  the  laser  beam,  and  hence  the  time 
available  for  the  ionization  processes. 

(U)   In  a  multispecies  environment  only  a  fraction  of  the  irradi- 
ance,  G,  say  a.G,  couples  to  each  species,  but  otherwise  the  analysis 
is  the  same,  and  the  emission  speed  scales  with  the  squareroot  of  the 
latent  heat  of  vaporization.  Actually  the  problem  is  even  more  compli- 
cated because  surface  species  will  be  heated  first;  so  the  steady  state 
emission  process  envisioned  here  is  an  idealization. 

(iD   Ignition  mechanisms  of  this  type  must  be  strongly  target 
material  sensitive.  They  will  also  be  sensitive  to  the  target  surface 
condition;  so  the  ignition  characteristics  might  change  from  pulse  to 
pulse  as  the  laser  flux  cleaned  the  surface.  This  pulse  to  pulse  vari- 
ation is  plausible  whether  the  actual  mechanism  depend  upon  an  oxide 
coating,  inclusions,  or  dirt.   If  the  laser  pulse  rate  is  not  too  high 
an  oxidized  surface  maintained  by  the  atmosphere,  again  Al  is  a  good 
example,  would  have  time  to  reconstitute  between  pulses. 

(U)   Because  this  mechanism  depends  upon  pure  heating  it  would 
depend  strongly  upon  the  laser  intensity,  but  rather  weakly  upon  the 
laser  frequency. 
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(U)  Very  high  energy  neutral  particles  are  also  emitted  from 
targets7.  The  energy  of  CO  molecules  have  been  measured  by  time-of- 
flight  techniques  to  be  14  eV  at  10  MW/cm2  irradiance  and  500  eV  at 
140  MW/cm2  from  Q-switched  ruby  lasers.  Thermal  emission  at  these 
energies  would  require  a  surface  temperature  of  the  order  of  105  K, 
which  is  highly  unlikely,  because  the  CO  molecule  would  certainly  dis- 
sociate at  such  temperatures.   It  is  possible  that  the  neutral  was  formed 
by  recombination  outside  the  target  after  emission,  or  that  these  high 
energy  neutrals  were  formed  by  charge  exchange.  These  effects  have  been 
included  here  for  completeness,  but  there  is  only  a  low  probability  that 
these  high  energy  neutral  species  are  directly  emitted  from  the  target; 
so  they  are  not  candidates  for  ignition  source. 

2.12  (u)   Ion  emission  adds  an  additional  complication  to  the  atom  emis- 
sion discussed  above.   If  the  surface  temperature  is  sufficiently  high 
an  appreciable  fraction  of  the  emitted  atoms  will  be  ionized.   If  the 
emitted  particles  are  singly  ionized  the  Richardson-Smith  equation 
applies,  and  the'  ion  curre'nt  is  given  by: 

j+  =  A+T2  exp(-0+/kT) 

where  A  is  a  constant,  T  is  the  absolute  temperature,  and  <f>     is  the 
positive  ion  work  function.  <f>     is  essentially  unknown,  but  it  can  be 
estimated  by  the  argument  that  the  energy  required  to  remove  the  singly 
charged  ion,  <f>   ,  plus  the  energy  required  to  remove  a  single  electron, 
0,  should  equal  the  energy  required  to  remove  a  neutral  atom,  £  ,  plus 
the  ionization  energy  of  the  neutral  atom,  I  .  This  argument  applied 

EL 

to  Al  gives  fzS  (Al)  =  5.3  eV,  and  to  Cu  gives  0  (Cu)  =  7.25  eV.   For 
comparison,  0(A1)  =  3.6  eV  and  0(Cu)  =  4.4  eV.  Thus  the  positive  ion 
emission  will  be  considerably  smaller  than  the  electron  emission  at  the 
same  temperature.  The  Richardson-Smith  equation  is  only  an  approxima- 
tion; so  it  should  be  used  with  caution. 

(U)  Even  less  reliable  is  the  Saha-Langmuir  equation  which  is 
derived  under  equilibrium  assumption,  but  is  often  applied  to  determine 
the  fraction  of  the  particle  current  which  is  ionized  by  the  laser  beam. 
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In  particular,  target  surface  temperatures  obtained  in  this  way  are 
highly  suspect.  The  Saha-Langmuir  equation  predicts  that  the  fraction 
ionized  takes  the  form 

j+/Jo-eXp[(^-Ia)/kT]. 

As  the  ionization  potential  is  generally  larger  than  the  work  function 
the  degree  of  ionization  must  increase  rapidly  with  temperature.  From 
Q-switched  ruby  laser  data  at  6  MW/cm2  the  duration  of  the  emitted  ion 
pulse  is  ~  100  nsec,  which  is  of  the  order  of  magnitude  required  for 
the  prompt  ignition  of  LSD  waves . 

(U)  Time-of-flight  energy  measurements  on  the  ions  emitted  by 
Q-switched  ruby  laser  experiments  give  energies  of  the  order  of  1000  eV. 
As  surface  temperatures  never  approach  such  values  plasma  effects  must 
be  influencing  the  ion  emission  process.  Multiphoton  effects  of  this 
magnitude  are  considered  to  be  totally  out  of  the  question. 

(U)   Ion  emission  from  oxide  layers  and/or  'dirt'  is  certainly  pos- 
sible, but  is  too  dependent  upon  the  details  of  a  particular  targets 
condition  for  general  discussion.  High  energy  molecules  containing  atoms 
from  surface  species,  particularly  oxides,  are  readily  detected  as 
molecular  ions,  but  the  ionization  may  have  occurred  after  emission. 

(u)  The  preceding  discussion  makes  use  of  temperatures,  which  im- 
plies that  the  materials  are  at  least  approximately  in  thermal  equilib- 
rium. This  assumption  may  not  be  justified,  even  though  no  non-equilib- 
rium effects  have  ever  been  positively  identified.   In  particular,  in 
the  presence  of  the  electron  emission  processes  discussed  in  the  next 
section,  strong  electrostatic  forces  must  be  applied  to  ions  near  the 
target  surface,  and  non-equilibrium  emission  processes  become  possible. 
When  many  electrons  have  been  stripped  from  the  surface  region,  heavy 
positively  charged  ions  remain;  so  that  a  large  fraction  of  the  target 
atoms  emitted  could  be  ionized.  Charge  conservation  at  the  target  sur- 
face becomes  a  major  consideration. 

(U)   In  summary,  it  is  possible  to  visualize  processes  by  which  the 
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target  could  emit  a  significant  number  of  ions  in  a  short  period  of  time, 
but  the  processes  depend  upon  detailed  considerations  of  particular  mater: 
ials  and  surface  conditions  which  defy  analysis  at  the  present  time. 

2.13  (U)  Electron  emission  processes  are  the  most  likely  ignition  mech- 
anism for  LSD  waves.  Thermionic  emission  is  known  to  be  present,  and  if 
the  laser  pulse  shape  is  known  the  shape  of  the  electron  pulse  can  be 
calculated  from  the  Richardson-Dushman  equation.  This  equation  can  be 
derived  from  either  thermodynamic9  or  kinetic1^  theory  arguments.  The 
usual  form  is 

j_  =  A_T2  expC-0/kT), 

where  j  is  the  electron  density,  0  the  electron  work  function  of  the 
target  surface,  and  A  is  a  universal  constant  (A  =  60.2  A  cm  2deg  2) 
for  many  metals.  Ready7  shows  that  the  Richardson-Dushman  equation  is 
a  reasonable  description  of  the  electron  emission  from  normal  laser 
pulses,  and  that  thermionic  emission  of  the  order  of  mA/cm2  is  readily 
produced  within  50  nsec  by  Q-switched  ruby  lasers11.  These  measurements 
were  performed  at  10  8  torr.  At  higher  pressures  where  absorbed  gases 
become  more  prevalent  the  electron  emission  rises,  which  indicates  that 
electron  emission  can  also  occur  from  surface  materials. 

(U)  Thermionic  emission  is  also  influenced  by  surface  morphology. 
If  the  surface  contains  pits,  protrusions,  or  inclusions,  the  local 
temperature  may  be  much  higher  than  would  be  predicted  from  the  laser 
flux  and  thermal  properties  of  the  bulk  target  material.  The  tempera- 
ture at  the  tip  of  a  protrusion  will  be  significantly  higher  than  that 
of  the  surrounding  material  because  the  tip  is  unable  to  dissipate  heat 
to  the  bulk  of  the  material  so  rapidly  as  a  planar  surface  can. 

(U)   Surface  atoms ,  whether  oxide  layers,  or  'dirt',  disrupt  the 
lattice  periodicity  of  the  target,  and  affect  the  work  function  of  the 
surface.  At  0  =  5  eV,  a  one  percent  reduction  of  the  work  function 
causes  an  80  percent  increase  in  j  at  1000  °K.  A  few  surface  atoms 
produce  noticable  effects,  a  monolayer  of  coverage  completely  changes 
the  work  function.  Because  all  materials  emit  thermionically,  at  high 
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irradiance,  target  material  effects  will  appear  primarily  in  the  time 
dependence  of  j  ,  that  is  in  the  shape  of  the  electron  pulse  ignited  by 
the  laser  pulse. 

(U)  Field  enhanced  emission  from  cold  targets  is  a  very  prompt 
mechanism,  and  therefore  is  attractive  as  a  LSD  wave  ignition  process. 
Musal  has  explored  this  possibility12  from  the  viewpoint  of  the  Fowler 
and  Nordheim  analysis13  which  is  also  discussed  in  the  Dushman  article 
cited  above9.  Musal  has  done  the  calculations  for  laser  irradiances  in 
the  range  from  50  to  500  MW/cm2  and  finds  that  the  electric  field 
strength  is  approximately  two  orders  of  magnitude  too  small  for  the  effect 
to  be  significant  for  electron  emission  from  a  plane  surface.  However, 
when  surface  features  are  included  in  the  analysis  the  field  strength  is 
greatly  enhanced,  and  the  mechanism  becomes  feasible.  Photon  assisted 
field  emission  has  recently  been  discussed  by  Caroli,  et.al.14. 

(U)  Schottky  emission,  thermionic  emission  in  the  presence  of  an 
external  electric  field,  has  also  been  extensively  studied.  This  is  a 
quantum  mechanical  barrier  penetration  process  in  which  the  work  function 
is  field  dependent15.  The  analysis  shows  that  the  work  function  will  be 
reduced,  but  that  otherwise  the  process  is  similar  to  thermionic  emission15 
Fields  of  the  order  of  105  volts/cm  cause  a  very  small  lowering  of  0.  The 
power  of  T  in  the  pre-exponential  factor  of  the  Richardson-Dushman  equa- 
tion is  difficult  to  determine  experimentally,  particularly  in  the  field- 
enhanced  environment.  The  factor  may  well  be  T3  rather  than  Richardson's 
T2  in  the  LSA  experiments,  because  there  is  considerable  experimental 
simularity  between  the  ignition  experiments  and  the  'complete  photo- 
electric effect'16. 

(U)  Multiphoton  emission  processes  can  produce  significant  electron 
current  density  from  ruby  and  Nd  doped  glass  laser  photons.  The  mechanism 
is  hardly  feasible  from  the  low  energy  photons  of  the  CO2  laser.  The 
mechanism  has  been  included  here  because  it  gives  the  earliest  response 
of  any  of  the  mechanisms  proposed.  Multiphoton  processes  become  partic- 
ularly significant  for  pico-second  laser  pulses  when  the  irradiance  is 
small  so  that  there  is  little  temperature  rise  in  the  target.  Thermionic 
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emission  decreases  exponentially  as  the  irradiance  decreases,  while  a 
two-photon  emission  process  decreases  as  the  square  of  the  irradiance; 
so  at  some  level  the  two  effects  must  be  equal,  but  this  normally  occurs 
at  irradiance  levels  well  below  those  required  to  initiate  LSA  waves. 

2.14  (U)  Photon  emission  mechanisms  are  poor  candidates  for  LSA  ignition. 
The  two  immediate  possibilities  are  reflected  photons  from  the  laser 
pulse  and  photons  generated  in  the  target.  The  reflected  photons  are  not 
sufficiently  energetic  to  cause  appreciable  ionization,  but  their  effect 
can  be  included  under  the  field  emission  heading  because  they  do  increase 
the  effective  electric  field12.  Photons  from  the  target  could  be  either 
thermally  emitted  by  the  hot  target  surface,  or  emitted  as  the  result  of 

a  multiphoton  atomic  process  which  occurs  at  the  surface. 

(U)  Although  sufficient  energy  is  available  to  melt  the  target  sur- 
face and/or  vaporize  it  under  some  conditions,  the  temperatures  reached 
are  realtively  low,  and  only  low  energy  photons  would  be  produced.  Such 
photons  probably  exist,  but  they  can  have  little  effect  upon  the  transi- 
tion from  ignition  to  plasma;  so  it  seems  unlikely  that  they  can  be  major 
contributors  to  the  LSA  initiation  process. 

(U)  Multiphoton  atomic  process  could  make  some  contribution  for 
ruby  laser  or  Nd/glass  illumination,  but  not  for  CO2  illumination  because 
the  photon  energy  is  so  small  that  10-20  photon  process  would  be  re- 
quired. There  is  not  experimental  evidence  that  photon  processes  contrib- 
ute significantly  to  LSA  ignition. 

2.15  (U)  Thermomechanical  response  of  the  target  to  the  laser  pulse  has 
been  proposed  as  an  ignition  mechanism  by.  Steverding12.   In  this  model, 
rapid  heating  in  a  shallow  region  near  the  target  surface  produces  a 
compression  stress  wave  which  propagates  to  the  target  surface  as  a 
shock  wave.  The  shock  wave  produces  spallation  at  the  surface.  The 
emission  velocity  is  the  shock  velocity,  and  the  emitted  material  is  a 
plasma  because  the  shock  wave  has  both  vaporized  and  ionized  the  target 
material. 
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(U)  As  envisioned,  the  target  is  heated  for  a  relatively  long  time 
(20  ysec)  and  the  laser  energy  is  deposited  in  a  very  shallow  region 
(10_1+cm  =  1  urn)  .  This  last  assumption  requires  a  depth  of  penetration, 
a  skin  depth,  of  approximately  X/10  for  the  CO2  laser  pulse.  A  region 
of  this  depth  is  consistent  with  the  calculation  of  Ready11,  for  the 
first  50  nsec  of  the  pulse,  but  thermal  conduction  increases  the  depth 
by  an  order  of  magnitude  within  the  first  200  nsec.  This  mechanism 
appears  to  be  too  slow  for  LSD  ignition,  but  it  may  be  feasible  for  LSC 
ignition. 

2.2  IGNITION  EXPERIMENTS 

(U)  Experimental  study  of  LSA  wave  ignition  processes  cannot  pres- 
ently be  separated  from  effects  attributable  to  transition  processes, 
but  certain  mechanisms  can  be  ruled  out  by  more  general  properties  of 
LSA  waves.  Here  some  of  the  experimental  properties  of  LSC  and  LSD 
waves  are  summarized  and  interpreted  in  the  light  of  the  various  pro- 
posed mechanisms.   Ignition  of  LSC  waves  differs  sufficiently  from  LSD 
ignition  that  they  are  considered  separately.  Much  more  experimental 
effort  has  been  devoted  to  LSD  waves,  see  the  reports  by  Hall,  et.al.18; 
Walters6;  and  Musal12;  so  the  ignition  properties  are  more  completely 
determined.  Pertinent  LSC  results  have  been  reported  by  Fowler  and  Smith19; 
Conrad,  Mangum,  and  Gurley20;  and  Thomas21  who  has  used  computer  simula- 
tions to  good  effect. 

(U)   This  survey  will  not  discuss  LSB  (Laser  Supported  Blast) 
effects  in  detail.  These  spherical  shock  waves  occur  in  conjunction  with 
LSC  and  LSD  waves,  but  are  not  supported  by  the  laser.  Hall,  et.al.18 
have  pictures  and  a  discussion.  Target  vapor  'microjets'  are  also 
observed18,  which  move  normal  to  the  target  surface;  so  if  the  target 
is  inclined  with  respect  to  the  laser  beam  direction  they  are  easily 
distinguished  from  LSC  and  LSD  waves.  Microjets  'ignite'  much  more 
slowly  than  LSD  waves,  but  on  a  time  scale  comparable  to  LSC  wave 
ignition  times. 
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(U)  Vacuum  ignition  of  Al  targets  by  nsec  duration  Nd-glass  laser 
pulses  of  irradiance  G  ~  1010  W/cm2  have  been  studied  in  some  detail 
by  Brooks22.  Note  that  the  targets  were  held  in  a  relatively  good 
vacuum  (~  10  5  torr) ,  and  the  irradiance  is  three  orders  of  magnitude 
greater  than  that  associated  with  LSD  wave  ignition,  but  that  the  total 
energy  absorbed  by  the  target  is  comparable  to  that  reported  in  the  LSC 
and  LSD  experiments . 

2.21  (U)   Laser  supported  combustion  wave  ignition  occurs  in  the  irradi- 
ance range  2  x  105  <_ G  <  6  x  106  W/cm2.  Hall,  et.al.18  has  observed 
LSC  waves  in  the  short  pulse  (~  20  ysec)  C02  laser  system  in  which  their 
LSD  experiments  were  performed.  They  have  not  made  extensive  studies 
of  LSC  wave  initiation  or  propagation.  The  detailed  information  currently 
available  comes  from  the  cw  CO2  laser  studies  of  Fowler  and  Smith19,  and 
the  long  pulse  (~  1  sec)  CO2  laser  studies  of  Conrad,  Mangum,  and  Gurley20, 

(U)  When  compared  with  the  LSD  wave  experiments  the  single  most 
outstanding  difference  is  the  relatively  long  initiation  time  of  the  LSC 
waves.   If  the  initiation  time  is  defined  as  the  time  until  the  LSA 
becomes  detectable,  both  ignition  and  transition  must  occur  during  the 
initiation  period.  This  time  is  of  the  order  of  200  ysec  for  LSC  waves, 
which  must  be  compared  with  LSD  initiation  times  of  the  order  of  50  nsec. 
LSC  is  very  much  slower,  more  than  three  orders  of  magnitude,  so  the 
ignition  processes  must  be  very  different.  On  the  assumption  of  a  linear 
rise  and  that  G  occurs  after  200  ysec,  the  target  surface  absorbs 

,     P  CO 

~  100  J/cm2  from  a  G  =  106  W/cm2  laser  before  initiation.  With  this 

P 
time  span  and  this  much  energy  available,  even  the  slower  ignition  mech- 
anism, such  as  atom  and  ion  emission,  are  feasible.  When  comparable 
laser  irradiations  are  applied  to  targets  in  vacuum  chambers,  where  the 
pressure  is  of  the  order  of  0.1  atm,  plasma  micro jet  phenomena  are 
observed.  Thus  target  vaporization  with  consequent  emission  of  atoms 
and/or  ions  seems  to  be  an  observable  characteristic  of  this  irradiance 
region.  Thomas21  has- done  computer  simulations  based  upon  vaporization 
ignition  which  compare  very  favorable  with  the  MSNW  experiments23. 
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(u)   In  summary,  target  vaporization  mechanisms  appear  to  be  suffic- 
iently to  explain  the  existing  information  about  LSC  ignition,  but  the 
necessity  of  these  mechanisms  has  not  been  established.  The  faster 
electron  emission  mechanisms  certainly  have  the  opportunity  to  function, 
but  the  irradiance  may  be  insufficient  to  initiate  a  surface  plasma. 
Sufficient  energy  has  been  dumped  into  the  surface  that  thermionic 
emission  must  occur  but  the  field  strength  contribution  from  the  irradi- 
ance minimized  the  Schottky  emission  contribution.  Thus  there  is  indirect 
evidence  that  the  Schottky  mechanism  is  an  important  contributor  to  LSD 
ignition. 

2.22  (U)   Laser  supported  detonation  wave  ignition  occurs  in  the  irradi- 
ance range  2  x  107  <_  G  <_  5  x  107  W/cm2.  Most  of  the  pertinent  experi- 
mental information  on  LSD  ignition  has  been  obtained  by  Hall,  et.al.*8, 
and  by  Walters6  who  has  also  done  a  detailed  integration  of  the  experi- 
mental evidence  applicable  to  the  ignition  problem  available  at  the  time 
of  his  review. 

(U)   LSD  waves  ignite  very  promptly;  the  initation  time  is  of  the 
order  of  50  nsec.   Quite  independent  of  the  target  material  the  initiation 
occurs  after  the  target  has  absorbed  approximately  1.7  J/cm2  of  energy. 
This  value  was  obtained  by  Walters6  who  integrated  his  experimental  laser 
pulse  shape  out  to  the  initiation  time.  This  is  insufficient  energy  for 
uniform  surface  vaporization  to  occur,  but  it  is  consistent  with  either 
a  thermionic  or  Schottky  emission  mechanism.  At  this  irradiance  the  rms 
electric  field  strength  is  insufficient  to  produce  field  emission  from  a 
plane  surface;  so  the  Schottky  mechanism  is  indicated. 

(U)  Walter's  careful  studies  of  the  ignition  process  indicate  that 
LSD  waves  ignite  from  'numerous  individual  luminosity  sites'.  The  wave 
front  forms  when  the  emission  from  these  sites  combine  in  front  of  the 
target  surface.  He  has  done  careful  scanning  electron  microscope  (SEM) 
studies  of  target  surface,  and  has  compared  the  appearance  of  an  indi- 
vidual surface  feature  before  and  after  LSD  wave  ignition.  The  pictures 
are  in  ref.  6.  A  feature  identified  before  the  laser  pulse  is  immediately 
recognizable  after  the  pulse.  The  SEM  scans  an  area  50  by  250  ym  at  a 
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2000X  magnification.  Although  absolute  correlation  between  surface 
features  and  luminosity  sites  has  not  been  accomplished  the  evidence 
that  features  cause  luminosity  sites  is  very  strong. 

(U)  The  laser  pulse  does  cause  surface  damage,  and  it  does  change 
the  surface  features,  but  there  is  no  gross  change  in  any  feature  which 
could  be  associated  with  significant  vaporization  of  its  surface.  When 
all  of  the  data  are  taken  together  electron  emission  from  surface  fea- 
tures appears  to  be  the  most  likely  ignition  mechanism.  The  ignition 
may  be  thermionic,  field  emission  as  proposed  and  discussed  by  Musal12, 
or  a  Schottky  effect  combination  of  the  two  at  the  surface  feature. 
Vaporization  from  surface  features  can  certainly  occur,  but  it  must  be 
a  second  order  effect  because  the  features  exhibit  so  little  damage, 
and  electron  emission  is  dominant. 

(U)  Tables  I  and  II  from  Hall,  et.al.18  give  LSD  ignition  thresholds 
as  measured  by  image  convertor  photographs  and  by  measurement  of  a  re- 
flected signal.  The  agreement  is  relatively  good  when  the  lack  of 
uniformity  of  the  surfaces  is  considered. 

(U)  Table  III  is  Walter's  summary6  of  the  relative  importance  of 
the  various  LSD  ignition  models  which  have  been  proposed. 
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TABLE  I 

THRESHOLDS  FOR  IGNITION  OF  LASER-SUPPORTED  ABSORPTION  WAVES  USING 
A  28  CM  FOCAL  LENGTH  LENS  AND  USING  A  NEW  SPECIMEN  FOR  EACH  TEST 
AT  ONE  ATMOSPHERE.   THRESHOLD  DETERMINED  BY  IMAGE  CONVERTER  PHOTO- 
GRAPHS OF  LUMINOSITY. 


Material 

Copper  Foil 

Lucite 

Polished  Tungsten 

Copper  Plate 

Tungsten  Plate 
(unpolished) 

Aluminum 
(7075  alloy) 

Aluminum  Foil 

Tantalum 

Carbon 

Aluminum  Black 

(very  heavily  anodized  Al) 

Lead  Foil 
(polished) 

Lead  Plate 
(unpolished) 


E_  (joules/cm2)    q    (watts/cm2) 


266 

3.83 

X  107 

258 

3.72 

X  107 

229 

3.30 

X  107 

192 

2.8 

X  107 

186 

2.68 

X  107 

176 

2.53 

X  107 

160 

2.30 

X  107 

128 

1.84 

X  107 

88 

1.27 

X  107 

87 

1.25 

X  107 

64 

.92 

X  107 

43 

.62 

X  107 
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TABLE  II 
THRESHOLDS  FOR  IGNITION  OF  LASER- SUPPORTED  ABSORPTION  WAVES  USING 
A  28  CM  FOCAL  LENGTH  LENS  AND  USING  A  NEW  SPECIMEN  FOR  EACH  TEST 
AT  ONE  ATMOSPHERE  THRESHOLD  DETERMINED  BY  MEASURING  REFLECTED 
SIGNAL. 


Material 
Fused  Silica 
Lucite 
Lexan  (polycarbonate) 


H20 


Teflon 

Copper 

Painted  Aluminum 

(.001  cm  gloss  black  lacquer) 

Titanium 

Fiberglass  Epoxy 

Cork 

Nickel 

Tungsten 

Soda  Glass 

Painted  Aluminum 
(.004  cm  red  enamel) 

Aluminum  Foil 

Aluminum  7075  Plate 

Lead 

Alumina 

Painted  Aluminum 

(.004  cm  flat  black  lacquer) 


ET  (joules/cm2) 


q   (watts/cm2) 
nmax 


310 

4.5  X 

107 

310 

4.5  X 

107 

280 

4.0  X 

107 

220 

3.2  X 

107 

195 

2.8  X 

107 

191 

2.75X 

107 

173 

2.5  X 

107 

160 

2.3  X 

107 

158 

2.3  X 

107 

144 

2.1  X 

107 

137 

1.97  X  107 

124 

1.79  X  107 

106 

1.52  X  107 

85 

1.22  X  107 

84 

1.2  X 

107 

75 

1.08  X  107 

68 

.98  X  107 

59 

.85  X  107 

42 

.60  X  107 
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TABLE  III 

ASSESSMENT  OF  THE  ROLE  OF  VARIOUS  MECHANISMS  IN  LSD  WAVE 

INITIATION  ON  PRACTICAL  ALUMINUM  SURFACES 


Mechanism 


Role 


Experimental  Basis 


(1)  Uniform  target        None 
vaporization 

(2)  Oxide  vapor  Minor 
absorption 

(3)  Chemical  reaction  Unknown* 

(4)  Thermionic  emission    Major 


(5)  Planar  reflection 

No  inde- 

pendent 

role  for 

G  <5X108 
P 

(6)  Local  target  heating 

Major 

(7)  Contaminant  and  gas 

Variable 

desorption 

(8)  Field  emission 

Minor  for 

GP<3X108 

(9)  Shock  heating 

Unknown* 

of  the  air 

(lO)Nonequilibrium 

Unknown* 

ionization  of 

vapor 

(11) Defect -enhanced 

Minor  for 

fields 

G  <3X108 
P 

SEM,  photography,  breakdown 
time 

Oxide-free  surface  initiation 


Target  emission,  breakdown 
time,  pressure  dependence 
of  threshold 

Clean-air  breakdown 


SEM,  photography 
Contaminant -free  surface 


SEM,  bent  sample,  breakdown 
time 


Breakdown  time 


Spectroscopy  now  underway  should  provide  assessment 
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2.23  (u)  Plasma  emission  experiments  are  discussed  in  this  section. 
Although  there  is  no  direct  evidence  that  these  results  bear  directly 
upon  the  LSA  ignition  problem  the  similarities  are  sufficiently  striking 
that  the  connection  may  ultimately  be  made. 

(U)  The  experiments,  performed  by  Brooks22,  used  a  Q  switched  Nd- 
glass  laser  with  energy  between  2.5  and  12  joules.  The  pulse  duration 
was  25  nsec,  and  the  focal  spot  size  was  ^  0.018  cm2;  so  the  irradiance 
was  1.4  x  109  <_  G  <_  5.7  x  109  W/cm2.  Note  that  the  total  energy  fluence 
absorbed  by  the  target  (140-700  J/cm2)  is  comparable  to  that  reported  in 
the  LSC  and  LSD  experiments. 

(U)  The  laser  pulse  produces  a  burst  of  plasma  from  an  Al  target. 
Careful  analysis  of  the  plasma  by  probes  determined  that  it  consists  of 
3  components  which  are  described  in  Table  IV. 

TABLE  IV 

Emission 
Component Time  Speed Density 

first  early         3  nsec     1.1  x  106  m/sec    2.5  x  1012/cm3 

second  early       12  nsec     5.9  x  105  m/sec    2.5  x  1012/cm3 

main  0-10  nsec     1.1  x  105  m/sec         10llt/cm3 

(u)  A  very  prompt  (<  1  nsec)  electromagnetic  pulse  was  also  detec- 
ted. Experimental  studies  confirmed  that  its  wavelength  is  less  than 
0.4  y;  so  it  was  not  a  reflection  of  the  incident  laser  pulse. 

(U)  An  applied  magnetic  field  has  little  effect  on  the  early  pulses; 
so  they  cannot  be  purely  electronic  in  nature.  Heavy  ions  must  be  present; 
hence  the  term  "plasma  pulse"  seems  appropriate.  The  main  pulse  is  delayed 
by  the  applied  magnetic  field  (B  a  950  gauss) .  The  delay  is  consistent 
with  an  interpretation  which  finds  the  plasma  front  relatively  rich  in 
electrons  while  the  main  body  is  richer  in  ions24  2e. 

(U)  The  early  plasma  pulse (s)  have  been  mapped  with  probes.  Figure 
1  shows  the  plasma  extent  after  20  nsec.  The  two  early  pulses  have  not 
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separated  at  this  time.  Note  that  the  plasma  is  asymmetrical  and  appears 
to  be  expanding  along  the  specular  reflection  direction.  The  plasma  has 
expanded  out  to  the  beam  spot,  and  is  now  eroding  the  entire  target  sur- 
face. Experiments  performed  on  large  area  targets  show  that  after  50 
shots  the  naked  eye  visual  erosion  region  has  an  area  of  —  I  cm2  when  the 
focal  area  is  0.018  cm2.  Microscopic  analysis  of  a  larger  region  has  not 
been  carried  out,  Figure  1  suggests  that  the  plasma  extends  beyond  this 
range . 

(u)  The  plasma  contours  after  160  nsec  are  shown  in  Figure  2. 
Note  that  the  two  early  pulses  have  now  separated  from  the  main  plasma 
(one  is  still  visible  on  the  right  hand  side  of  the  figure) .  The  plasma 
emission  is  now  symmetrical,  and  the  main  pulse  is  moving  normal  to  the 
target  surface. 

(u)  The  situation  after  800  nsec  is  shown  in  Figure  3.  Note  in 
particular  that  the  plasma  is  still  attached  to  the  target;  so  that  it 
is  still  feeding  upon  the  target  by  sputtering. 

(u)  This  material  has  been  included  because  it  demonstrates  the 
target  material  behavior  when  it  expands  into  an  atmosphere  of  comparable 
density  (n  CT  1013/cm3).  At  a  standard  atmosphere  the  plasma  would  inter- 
act much  more  strongly  with  the  background  gas;  so  the  development  must 
be  much  more  complicated  than  illustrated  here.  Still  certain  similari-'  ■ 
ties  should  remain.  The  prompt  emission  from  the  target  should  not  be 
background  pressure  sensitive.  The  electron-rich  front  surface  of  the 
main  plasma  burst  is  a  target  phenomenon  which  should  carry  over.  Note 
that  it  supports  the  idea  that  the  initial  emission  must  be  electronic 
with  ions  following  closely.  This  effect  will  be  discussed  in  greater 
detail  below.  Even  though  the  laser  pulse  has  been  extinguished  for 
almost  0.8  ysec  in  Figure  3  the  plasma  is  still  feeding  on  the  target 
as  it  expands  away  from  the  surface.  There  is  strong  evidence  that  the 
radial  expansion  of  the  plasma  is  constrained  by  a  self-induced  magnetic 
field22"27.  This  would  help  to  explain  the  long  time  of  contact  between 
the  plasma  and  the  target  surface. 
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Fig.  2  Relative  density  contours  for  a  laser  produced 
plasma  160  nsec  after  ignition. 
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produced  plasma  800  nsec  after  ignition. 
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(U)  Electron  microscope  examination  of  the  focal  area  on  the  tar- 
get produces  evidence  for  shock  wave  damage  of  the  target.   "Large" 
masses  of  Al  have  been  spalled  from  the  target.  There  is  also  evidence 
for  "splashing"  of  liquid  Al  from  the  target  onto  the  probes  when  they 
are  located  close  to  the  target  surface.  At  these  very  high  irradiances 
and  with  these  very  short  pulses,  Steverding's  model17  seems  to  apply 
during  some  part  of  the  ignition  event. 

(y)  These  investigations  also  have  a  bearing  upon  the  nature  of 
the  electron  emission  mechanism  which  causes  prompt  ignition.  The 
electron  energies  corresponding  to  the  emission  speeds  of  the  plasma 
components,  see  Table  IV,  are  first  early  =3.4  eV,  second  early  =  1  eV, 
main  =  3  x  10  2  eV.  These  results  are  consistent  with  an  interpretation 
which  suggests  that  the  first  two  components  are  caused  by  a  high  field 
emission  mechanism,  while  the  main  component  is  thermionic.  This  inter- 
pretation has  been  discussed  by  Thomas,  Musal  and  Chou28.  The  structure 
of  the  components  suggests  that  the  electrons  are  emitted  first  and  that 
ions  from  the  target  surface  are  then  pulled  along  behind  the  electrons 
by  electrostatic  forces. 
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3.0  QUALITATIVE  SUMMARY  OF  TRANSITION  PROCESSES 

(U)  The  ignition  process  is  complete  as  soon  as  conditions  have  been 
established  in  the  air  in  front  of  the  target  surface  which  will  ultimately- 
lead  to  the  production  of  an  LSA  wave.  A  certain  number  of  charged  par- 
ticles must  be  present  (both  electrons  and  ions)  and  the  number  of  both 
must  increase  with  time.  The  transition  to  full  LSA  behavior  is  extremely 
complex  because  a  large  number  of  processes,  both  productive  and  dissi- 
pative,  are  occurring  simultaneously.   If  an  LSA  wave  is  to  result,  the 
productive  mechanisms  must  dominate.  Any  mechanism  which  increases  the 
absolute  number  or  density  of  electrons  is  productive,  any  which  decrea- 
ses the  electron  population  is  dissipative.  Attention  can  be  focused  on 
the  electrons  because  an  electron  produced  is  also  an  ion  produced.  The 
converse  of  this  statement  is  not  necessarily  true  because  electrons  can 
diffuse  out  of  the  ionized  region  when  it  is  small.   In  most  cases  this 
surface  leakage  is  not  a  dominant  consideration,  but  it  is  a  significant 
contributor  when  the  dimensions  of  the  ionized  volume  approach  the  mean 
free  path  of  electrons  in  air.  This  number  depends  upon  many  considera- 
tions, but  is  of  the  order  of  0.5  ym  in  STP  air.   In  general  focal  spots, 
and  hence  ignition  plasmas  are  much  larger;  so  electron  diffusion  is  a 
secondary  consideration. 

(U)  The  atomic-scale  production  and  dissipation  processes  are  sum- 
marized in  Figure  4.   In  the  following  material  A  will  stand  for  an 
unexcited  atom  or  molecule,  A*  for  an  atom  in  an  excited  state  and  A+ 
for  the  singly  charged  ion  of  that  atom.  Multiple  ions,  A2+  etc.,  and 
the  transition  from  the  molecular  to  the  atomic  state  will  be  ignored  in 
this  discussion.   Both  complicate  the  already  complex  discussion  without 
significantly  contributing  to  the  analysis.  The  atom-molecule  transition 
is  an  additional  loss  mechanism  which  reduces  the  effectiveness  of  the 
production  mechanisms  discussed  below. 

(U)  All  of  the  productive  processes  feed  directly  or  indirectly  upon 
the  energy  which  the  laser  pulse  supplies  to  the  ionized  region.  This 
region  will  be  referred  to  as  the  plasma,  even  though  in  the  strictest 
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Fig.  4  Schematic  representation  of  production  and 
dissipation  processes.  Arrows  to  the  right  represent 
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sense  absolute  charge  neutrality  can  not  be  guaranteed  because  of  electron 
diffusion  and  target  surface  effects. 

(U)   Because  photons  couple  poorly  to  atomic  mass  particles  (the  mass 
mismatch  is  too  large)  photon-electron  processes  play  the  primary  role. 
These  processes  can  be  catagorized  as  bound-bound,  bound-free,  and  free- 
free  where  the  descriptive  words  apply  to  the  condition  of  the  electron 
before,  and  then  after  the  photon-electron  collision.  Although  these 
processes  must  occur  'near'  an  atom  or  ion  to  conserve  momentum,  the 
photon- electron  collision  mechanism  dominates  the  event.  The  processes 
are 

1)  hv  +  A  ->  A*  (bound-bound) 

2)  hv  +  A  (A*)  ■*■  A  +  e~  (bound- free) 

3)  hv  +  e  -*■  e  (free-free) 

Mechanisms  1)  and  2)  combine  to  produce  multiphoton  ionization,  or  on  a 

macroscopic  scale,  microwave  heating  of  the  plasma28.  The  two  methods  of 
analysis  are  formally  equivalent.  This  is  a  plausible  plasma  heating 
mechanism  where  the  laser  photon  energy  exceeds  1.0  eV  (Ruby,  Nd-Glass, 
Nd-YAG)  but  it  is  not  practical  for  the  CO2  laser  photon,  which  is 
approximately  0.1  eV. 

(U)   In  the  absence  of  multiphoton  ionization,  mechanism  3)  must  be 
the  dominant  coupling  between  the  laser  pulse  and  the  plasma. 

(U)  Given  a  supply  of  energetic  electron  from  mechanism  3),  addit- 
ional ionization  can  be  produced  by  electron-atom  mechanisms.  The 
productive  mechanism  is: 

4)  e"  +  A(A*)  -*■  A  +  2e~.  (ionization) 

Note  that  ionization  will  be  facilitated  by 

_  1 

5)  e~  +  A  ->  A*  +  e  .  (excitation) 

Mechanism  4)  is  the  ultimate  plasma  source,  while  in  one  sense  mechanism 
5)  is  a  loss  mechanism  because  it  withdraws  energy  from  the  electron 
heating  process.   It's  exact  catagorization  will  depend  upon  the  details 
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of  a  particular  plasma  condition  and  upon  the  electron  and  ion  tempera- 
tures in  that  plasma. 

(u)  Two  additional  mechanisms: 

_  i 

6)  e  +  A*  ■+■  A  +  e  (collision  de-excitation) 


and 


7)  e  +  A  ■*■  A  +  hv'  (recombination) 


are  clearly  dissipative. 

(U)  A  great  deal  of  study  has  been  done  on  photon-electron  and 
electron-heavy  particle  scattering.  Much  of  the  effort  required  to  trans- 
late these  investigations  into  forms  useful  for  the  LSA  transition  analysis 
has  been  done  by  Ready7  (see  in  particular  the  sections  on  spark  ignition 
in  air)  and  by  Zel 'dovich  and  Raizer29  (particularly  in  Vol.  I,  which 
contains  the  fundamental  physical  material). 

( U)   It  should  be  noted  that  a  description  of  the  transition  from  the 
atomic  interaction  mechanisms  to  plasma  ignition  is  dependent  upon  some 
form  of  transport  theory,  which  in  turn  requires  local  (psuedo)  thermo- 
dynamic equilibrium.  The  electron  equilibration  time  is  such  that  a  well 
defined  electron  temperature  seems  plausible,  but  the  concept  of  an  ion 
temperature  should  be  used  with  great  caution,  because  much  of  the  ions' 
energy  may  be  associated  with  the  streaming  of  the  plasma  away  from  the 
target.  Also,  the  ion  equilibration  time  is  approximately  four  orders 
of  magnitude  greater  than  the  electron  equilibration  time,  which  places 
it  uncomfortably  close  to  the  entire  initiation  time  of  the  LSA  wave. 
In  summary,  the  electrons  probably  have  a  temperature,  the  ions  may  have 
a  temperature,  but  the  electrons  and  ions  are  clearly  not  in  thermal 
equilibirum. 

(U)  To  a  good  approximation  the  plasma  density  will  grow  exponentially 
with  the  number  of  high  energy  electron  "generations"  n,  required  to 
reach  a  given  level  of  ionization.  Thus 


N  =  N  exp 
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where  T  is  the  lifetime  of  one  electron  heating  cycle: 

_    _  i 

hv  +  e  •*■   e 


followed  by 


+   _ ' 
e  +  A(A*)  -v  A  +  e  . 


For  an  ignition  plasma  of  the  type  studied  by  Brooks22  approximately  16 
generations  are  required  to  produce  plasma  densities  approaching  100% 


ionization.   If  N  were  10  electrons,  43  generations  would  be  required, 
99%  of  the  ionization  occurs  in  the  last  7  generations.  The  breakdown 
process  can  proceed  in  times  as  short  as  a  few  nseconds29. 
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4.0  PROPAGATION  OF  LSA  WAVES 

(u)  The  preceding  two  sections  discuss  the  ignition  of  a  plasma  at 
the  target  surface  and  the  transition  of  that  plasma  to  a  high  density 
state.  This  section  is  devoted  to  the  behavior  of  the  system  once  the 
plasma  has  become  dense. 

(U)  The  transition  process  is  complete  once  the  plasma  reaches  a 
constant  density.   If  this  density  is  less  than  that  established  by  setting 
the  plasma  frequency  equal  to  the  laser  frequency 

N  e2 
e 

=  v  =  v. 


m      p    laser 


K 

v  =  8.97  x  10J  N   (cm  *) 


or 


N  (cm-3)  =  1.12  x  109/A2  (m) 

p v    '  p  *•  J 

the  plasma  is  transparent  to  the  radiation,  and  the  target  surface  continues 
to  supply  material  to  the  plasma.  Once  the  plasma  density  exceeds  the 
threshold  value  the  plasma  becomes  opaque  and  shields  the  target  from  the 
laser  pulse. 

(U)  Three  types  of  behavior  have  been  identified  experimentally.   When 
the  irradiance,  G,  is  of  the  order  of  104  W/cm2  a  plasma  is  formed  at  the 
target  surface,  but  its  density  never  exceeds  the  plasma  frequency  threshold, 

The  plasma  is  in  a  steady  state  under  the  combined  influences  of  the  laser 
irradiation,  convection,  reradiation  and  target  surface  erosion.  The 
resultant  PLASMATRON  is  transparent  to  the  laser  radiation;  so  it  remains 
attached  to  the  target  surface.   When  the  irradiance  is  of  the  order  of 
105  W/cm2  the  plasma  detaches  from  the  target  surface  and  propagates  rather 
slowly  (v  -   103  cm/s)  back  up  the  laser  pulse.  The  plasma  is  opaque  and 
protects  the  target  surface  from  the  laser  beam.  The  ignition  mechanism 
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in  the  air  in  front  of  the  plasma  is  very  complex.  This  low  velocity  dis- 
turbance is  the  LSC  wave.   Lastly,  when  the  irradiance  is  of  the  order  of 
107  W/cm2,  the  plasma  propagates  with  a  velocity  which  exceeds  the  speed 
of  sound  and  an  LSD  wave  results.  The  names  of  the  waves  come  from  their 
similarity  to  high  temperature  hydrodynamic  waves.  The  same  mathematical 
models  apply;  so  hydrodynamic  theory  can  be  used  to  describe  LSC  and  LSD 
waves.  This  analysis  was  first  developed  in  the  Soviet  Union  by  Raizer3 
and  in  the  West  by  Ramsden  and  Savic30. 

(U)  Detonation  waves  are  easier  to  discuss  theoretically  because  a 
number  of  simplifying  assumptions  can  be  made  in  the  analysis.  Combustion 
waves  are  more  complex  and  the  details  of  their  propagation  are  less  well 
understood.  The  analysis  here  passes  from  a  simple  one-dimensional  dis- 
cussion of  the  relevent  hydrodynamic  equations  to  a  presentation  of  the 
salient  characteristics  of  the  LSD  wave,  and  then  considers  the  LSC  wave 
where  closed  form  analytic  expressions  are  less  useful. 

4.1  HYDRODYNAMIC  EQUATIONS 

(U)  Consider  a  one-dimensional  system  of  the  form  shown  in  Figure  5, 
which  is  characteristic  of  the  LSD  wave.   Here  the  absorption  region  is 
optically  thick;  so  the  laser  radiation  is  absorbed  there.  The  wave  prop- 
agates to  the  right  with  speed  D.  The  region  denoted  "ambient  air",  will 
be  Region  1  in  what  follows,  and  is  "ahead  of"  the  wave.  The  region 
denoted  "relaxation  region"  will  be  Region  2,  and  is  "behind"  the  wave. 

(U)  The  analysis  presented  here  is  not  sophisticated,  but  it  includes 
all  of  the  physical  processes  and  variables  which  are  pertinent  to  this 
discussion.   It  differs  from  the  similar  hydrodynamic  derivations  because 
an  energy  source  term  has  been  included,  but  the  momentum  contribution 
from  the  laser  pulse  have  been  neglected.  A  more  detailed  discussion  is 
available  in  Zel'dovich  and  Raizer31.   In  the  material  which  follows  the 
thickness  of  the  discontinuity  (£  in  Figure  5)  will  be  neglected. 

(U)  The  variables  required  to  characterize  the  system  are: 
p  the  pressure 
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laser-supported-detonation  wave. 
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p  the  density  or  its  reciprocal 

V  =  1/p  the  specific  volume 

e  the  specific  internal  energy 

u  the  fluid  flow  velocity  in  the  LAB  system 

Occasionally,  h  =  e  +  PV  the  specific  enthalpy  will  be  used  to  simplify 
equations.   Here. the  term  "specific"  means  "per  unit  mass".  The  general 
term  "fluid"  has  been  used  rather  than  "air"  or  "gas"  because  the  material 
behind  the  discontinuity  is  a  plasma. 

(U)  The  analysis  considers  the  conservation  of  an  element  of  mass, 
Am,  its  momentum,  and  its  energy,  as  the  energy,  as  the  element  moves 
across  the  discontinuity.   In  terms  of  the  variables,  define  a  volume 
element,  6t,  on  each  side  of  the  discontinuity.  Then, 

6t  =  (D  -  ui)  •  A6t  =  (D  -  u2)  *  ASt; 

where  the  area  of  the  element  A,  remains  constant  as  it  passes  through  the 
discontinuity,  and  of  course  the  same  time  duration,  <5t,  applies, 

(U)  Then 

Am  =  pSx 

Am  =  pi(D  -  ui)  ♦  ASt  =  p2(D  -  u2)  •  ASt; 
which  gives 

Pl(D  -  ui)  =  p2(D  -  u2). 
This  conservation  of  mass  relation  is  most  conveniently  written  as 

D(P2  -  Pi)  =  P2u2  -  PlUl.  I' 

(u)   Conservation  of  linear  momentum  requires  that  the  momentum  change 

of  the  mass  element  be  equal  to  the  impulse  received.  Then 

-*       >      -> 
(P2  -  Pi)A<5t  =  Amu2  -  Amuj. 

-*   -A   _k  _i 

When  D,  A,  uj,  and  u2  are  all  parallel  this  equation  simplifies  to 

D(p2u2  -  p^i)  =  (P2  +  p2u2)  -  (Pi  +  Piu2).       II1 

2  1 

This  is  the  simplification  which  reduces  the  relation  to  a  one-dimensional 
form. 
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(U)  Conservation  of  energy  follows  from  the  First  Law  of  thermo- 
dynamics : 

AU  =  AQ  +  AW. 

The  element  of  mass,  and  hence  the  quantity,  U,  contains  energy  in  two 
forms:  the  specific  internal  energy,  e,  and  the  specific  kinetic  energy 
(l/2)u2.   (Watch  the  units!).  Then  one  writes 

AU  =  Am(e2  -  el) 
where 

e  =  e  +  (l/2)u2. 

The  heat  added  is  the  integrated  laser  flux  through  A  in  time  6t, 

AQ  =  G  A6t 
x    o 

and  the  work  done  on  the  element  of  mass  is 

AW  =  (P2u2  -  Piui)  -  A6t. 

When  all  terms  are  combined  the  relation  reduces  to 

D(P2£2-Pi£i)  =  G  +  (p2u2£2  +  P2U2)  -  (Piuiei  +  Piui).         Ill' 

(U)  These  equations  can  be  considerably  simplified  by  transforming 
to  a  coordinate  system  which  moves  with  the  discontinuity: 

u*  =  u  -  D. 

Note  that  e  /  e*  because  Of  the  kinetic  energy  term.  The  three  conser- 
vation equations  reduce  to: 

p:u*  =  p2u*  s  j, 

Pi  +  Piu*2  =  P2  +  P2U*: 
1  2 


and  U/2)u< 


(£2+P2V2+    -s-2-  )  -  (£i+PiVi+ 


I 

I 

II 

(l/2)u*2 

1      1    1 

G 
0 

=   0. 

TT 

2          }    ' 

P2U* 
2 
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It  is  convenient  to  define 
q  =  G  /p2u* 

O     2 

where  q  is  the  specific  energy  which  the  element  of  mass  has  acquired 
from  the  laser  pulse  as  it  moves  through  the  absorption  region,  see 
Figure  5. 

(U)  So  far  only  three  relations  have  been  specified  relating  four 
unknowns,  but  a  number  of  useful  results  can  be  obtained  before  the 
equation  of  state  of  a  specific  material  is  introduced.   The  equations 
are  most  easily  manipulated  when  written  in  terms  of  V.  With  this  modi- 
fication, Eq.  I  becomes 

u*/u*  =  V1/V2, 
1   2 

and  Eq.s  I  and  II  can  be  manipulated  to  obtain 


MP 

M  Vi-V: 


j     =    U*/V2    =\Jv—T7-      =    U*/Vl  (1) 

2  V  V1-V2  1 

These  relations  are  of  special  use  when  u\   =  0  (ambient  air  at  rest);  so 

that  u*  =  -D.   With  this  additional  condition 
1 


|D!  =  vAki' 


(2) 


this  is  the  speed  of  the  discontinuity  in  the  LAB  coordinate  system. 


(u)   It  is  possible  to  show  that  ui  <  D  if  the  discontinuity  is  to 
gate  as 
positive;  so 


propagate  as  indicated,  so  u*  must  always  be  negative.   But  Vi  must  be 


j  S  V(P2-Pi)/(Vi-V2) 


must  be  intrinsically  negative.  This  factor  will  always  be  assumed  to  be 
negative. 

(U)  When  ui  =  0,  Eq.  I  can  be  written  (use  Eq.  2) 
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U2  =  V 


(P2-PO/CV1-V2)  (3) 


N.B.  that  U2  is  in  the  LAB  system.  u2  <  D  for  all  possible  discontin- 
uities (otherwise  the  material  behind  the  discontinuity  would  catch  up 

to  it!),  so  u*  must  also  be  negative. 
2 

(U)   It  is  convenient  to  remove  the  microscopic  variables  u*2  and 

u*  from  Eq.  III.   From  Eq.  I 

u*2  .  u*2  =  -(V1+V2HP2-P1),   •  (4) 

2     1 

and  Eq.  Ill  becomes 

£2-£i  +  l/2(V2-Vi)(Pi+P2)  =  q  (5) 

The  negative  sign  in  the  definition  of  q  was  chosen  so  that  Eq.  5  could 
be  interpreted  physically.   It  is  consistent  with  the  fact  the  u*  is 
intrinsically  negative;  so  q  as  defined  will  always  be  positive. 

(U)  To  proceed  further  one  must  assume  an  equation  of  state  for 
the  material.  The  ideal  gas  will  be  used  in  this  discussion  because  it 
is  sufficient  for  demonstrative  purposes,  even  though  it  can  hardly  be 
expected  to  describe  a  plasma  correctly.   In  terms  of  the  specific  heat 
ratio  (C  is  the  specific  heat  at  constant  pressure,  C  at  constant 
volume) 

Y  S  Cp/Cv  ,  (6) 

the  ideal  gas  law  can  be  written  as 

P  =  (y-l)pe.  (7) 

At  times  it  is  more  convenient  to  define  a  temperature 
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T  E  e/Cy  (8) 

(note  that  C  is  assumed  to  remain  constant) ;  then 

P  =  (Y-l)pCvT.  (9) 

A  useful  auxiliary  relation  is 

(  ijjj  )PV  =  e  +  PV  =  h  =  CpT  (10) 

(U)  These  basic  relations  will  now  be  applied  to  various  LSA  wave 
propagation  conditions. 

4.2  ONE -DIMENSIONAL  LSD  WAVES 

(U)  Detonation  waves  occur  when  the  speed  of  the  wave  exceeds  the 
speed  of  sound  in  the  ambient  air,  D  >  ci.  The  wave  propagation  is  not 
stable  unless  it  also  meets  the  condition  U2  £  C2  £  D,  where  C2  is  the 
speed  of  sound  behind  the  discontinuity.   With  these  conditions  the  dis- 
continuity propagates  as  a  shock  wave,  and  many  of  the  analytic  results 
obtained  from  shock  wave  studies  also  apply  to  LSD  wave  propagation. 

(u)  The  system  is  easiest  to  analyze  if  the  CHAPMAN -JOUGET  condition 
is  applied: 


u*  =  c2, 


which, when  combined  with  the  ideal  gas  medium,  gives 

u*2  =  c  2  =  Y2P2V2.  (11) 

2     2 

Zel'dovich  and  Raizer31  discuss  the  meaning  of  this  condition  in  con- 
siderable detail. 

(u)  When  Eq.  11  is  used  with  Eq.  II,  one  obtains 
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P1V1  +  u*2 


in  conjunction  with  Eq.  I,  it  gives 

u*2  +  P1V1  y2 

V2=  (-Hi KifiTT'    '  (13) 

and  in  Eq.    Ill  with  Eq.    10  it  reduces  to 

Y2  Yi  t  u*2 

(  YTT-f  )p2V2   "    (  y—TJ  )pivi    -  7  V2P2V2   -     4~  =  q.  (14) 

Note  that  the  change  in  medium  properties  behind  the  shock  front  has 
been  accounted  for  in  this  analysis  by  allowing  Yi  ^  Y2 .   If  one  now 
eliminates  P2  and  V2  from  Eq.  14,  the  resultant  equation  is  a  bi-quadratic 
in  u*  .  The  solution  to  this  equation,  with  the  replacement 

P1V1  +  (Yi  -  DC  Ti 


is 


u*  =  |l/l/2(Y2-l)[(Y2  +  l)q  +  (Yi+Y2)CViTi]     j 

+     |Vl/2(Y2+l)[(Y2-l)q  +  CY2-Yi)CViTi]    J    .        (15) 

This  equation  becomes  very  simple  in  the  strongly  illuminated  condition 
where  q  >  >  C  Ti  (the  specifid  internal  energy  of  the  ambient  gas  is 
negligible  compared  to  the  specific  energy  added  by  the  laser).   It 
becomes 

1  1  r       2     -> l*2 

ju*J  =  JDJ=  [z(Y2   -  l)qj  (16) 
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Now  restate  Eq.  12  to  provide  the  pressure  ratio, 


P2 

pT 


u*2  +    (Yi   -   DCViTi 
(Y2+Yi)(Yi    -   l)CVlTi    ' 


(17) 


and  Eq.s  13  and  17  provide  the  volume  ratio 


V2 

v7 


Y2 


Y2+I 


u*z  +  (Yi  -  1)C  Ti 

1  vi 


u* 


(18) 


If  one  now  uses  the  requirement  (Yi 


Y2 


V2 

vT 


1)C  Ti  <  <  u*  ,  Eq.  10  reduces  to 
J   vi       1    n 


(19) 


which  applies  in  the  strong  shock  limit  of  a  strongly  fed  wave  which  also 
satisifies  the  CHAPMAN- JOUGET  condition.   To  the  same  approximation 


2(Y2-D 


Pi      (Yi-D   CViTi 


YiD' 


(Y2  +  DC- 


Yi 


Y2+1 


W 


(20) 


where  M  is  the  Mach  Number  of  the  detonation  wave  motion. 

(U)   If  one  replaces  q  by  its  definition  in  Eq.  16,  one  obtains  a 
relation  between  the  wave  speed  D  and  the  irradiance  of  the  laser: 


D  =    [2(Y2   -   l)ViG0] 


i/3 


(21) 
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Note  that  this  result  depends  upon  the  ambient  density  in  front  of  the 
shock,  the  specific  heat  ratio  behind  the  shock,  and  the  irradiance. 
The  relation  between  D  and  G  has  been  well  tested  experimentally. 

4.3  ONE -DIMENSIONAL  LSC  WAVE 

(U)  Combustion  waves  occur  when  the  velocity  of  propagation  of  the 
disturbance  is  less  than  the  velocity  of  sound.  Velocities  for  LSC 
phenomena  range  from  D  =  0,  the  plasmatron,  to  the  order  of  a  few 
meters/second.  Here  the  D  =  0  case  is  specifically  excluded. 

(U)  The  LSC  disturbance  is  a  more  complex  structure  than  is  the 
LSD  wave  as  sketched  in  Figure  6.  This  figure,  from  the  work  of  Jackson 
and  Neilsen32,  is  the  result  of  a  computer  computation.   Its  qualitative 
features  agree  very  well  with  the  experimental  data  available  on  LSC  and 
LSD  waves. 

(U)  The  LSC  wave  structure  is  conveniently  divided  into  three 
regions:  A,  the  ambient  air  in  front  of  the  disturbance;  B,  the  transi- 
tion region,  which  has  been  warmed  and  which  contains  some  ionized  par- 
ticles, but  which  is  transparent  to  the  laser  photons;  and  C,  the  plasma, 
which  is  opaque  to  the  laser  photons.  The  electron  density  and  tempera- 
ture increase  continuously  with  distance  from  the  wave  front  so  that  at 
the  plasma  front  the  plasma  frequency  of  the  ionized  medium  equals  the 
frequency  of  the  incident  laser  light. 

(U)  The  LSC  wave  differs  from  the  LSD  wave  because  the  disturbance 
created  by  the  wave  front  is  not  sufficiently  intense  to  produce  the 
electron  density  required  to  support  the  plasma;  so  additional  ionization 
must  be  generated  by  the  plasma.   For  the  LSD  wave  the  two  fronts  coalesce 
and  the  analysis  is  much  simpler. 

(U)   Because  of  the  more  complex  structure,  the  simple  conservation 
law  analysis  of  the  LSD  system  contained  in  Eq.s  I,  11^  and  III  is  no 
longer  appropriate.   It  could  be  applied  between  regions  A  and  C,  but 
most  of  the  interesting  phenomena  which  control  the  wave  development 
occur  in  region  B.   In  region  B  the  parameters  are  functions  of  distance 
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Fig.  6A  One-dimensional  LSD  wave  driven  by  a  C0_  laser  flux 

of  1  x  107  watts/cm2  decoupling  of  high  temperature  region 
and  shock  wave  does  not  occur  V-—.  =  4.0  x  105  cm/sec. 
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from  the  wave  front;  so  a  full  divergence  theory  analysis  is  required. 
In  addition,  the  functional  variations  are  driven  by  the  atomic  level 
process  discussed  in  Section  2;  so  the  resultant  non-linear  differential 
equations  are  extremely  complicated. 

(u)  Here  it  seems  appropriate  to  discuss  the  model  and  the  pheno- 
menological  parameters,  but  a  presentation  of  the  actual  numerical 
methods  used  to  solve  the  equations  would  go  beyond  the  intent  of  this 
survey.   The  interested  reader  will  be  referred  to  the  current  literature 
at  the  appropriate  points  in  the  discussion.   In  many  respects  this 
presentation  draws  heavily  upon  Nielsen  and  Canavan5,  but  there  are 
significant  variations  from  their  notation.  The  original  analysis  of 
LSC  waves  performed  by  Raizer33,  has  strongly  influenced  all  subsequent 
studies. 

(U)   Because  the  LSC  wave  propagates  so  slowly,  to  a  very  good 
approximation  the  pressure  may  be  assumed  constant  throughout  the  dis- 
turbance and  the  conservation  of  momentum  condition  contributes  nothing 
to  the  analysis.   Similarly,  there  is  very  little  mass  flow;  so  effec- 
tively only  the  energy  conservation  condition  remains.  This  condition, 
written  in  differential  form  is 

Div(pvh  +  Pv)  =  Q. 
Here  h  is  the  enthalpy  and  Q  contains  all  of  the  energy  source  and  sink 
terms.  Three  terms  can  be  identified  in  Q: 

Q  =  $  +  ¥  +  £  (22) 

where  $  is  the  radiation  term,  ¥  the  heat  conduction  term,  and  C  is  the 
electron  transport  term.   Each  of  these  terms  operates  somewhat  differently 
in  the  three  regions.  Q  vanishes  in  region  A,  because  by  definition  the 
medium  is  undisturbed  there.  The  laser  supplies  energy  in  region  C, 
because  the  plasma  is  opaque,  so  $  is  a  source  term  which  can  be  written 

$  =  K  G  (x)  (23) 

o  o v  J 

where  K  is  the  absorption  coefficient,  and  G(x)  is  the  residual  laser 
flux  at  distance  x  from  the  wave  front  which  is  taken  as  the  origin.  The 
laser  intensity  in  turn  must  satisfy 
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g«r«  (24) 

subject  to  the  boundary  condition  GCo)  =  G   the  laser  irradiance. 

o, 

(U)  The  absorption  coefficient,  K,  is  a  strongly  varying  function 

of  temperature,  see  Figure  7.   Because  the  coefficient  rises  so  sharply, 

in  most  calculations  it  is   approximated  by  the  function 

K(T)  =0.0        T  <  12,000°K 

=  0.7  cm'1   .  T  >  12,000°K  (25) 

Other  investigations  prefer  to  set  K(T)  =  0.6  cm"1  when  T  >  1.0  eV32. 

(u)   In  region  C,  the  thermal  conduction  term,  T,  is  a  loss  term. 
It  supplies  much  of  the  energy  which  increases  the  temperature  of 
region  B.  The  thermal  conductivity  of  air,  3,  is  an  extremely  compli- 
cated function  of  temperature,  see  Figure  8.  A  convenient  approximation 
for  3  is18 

3  =  2  x  10"6  T  (watts/cm  °K) .  (26) 

The  heat  flux  is  given  by 

d2T 
f  =  -B  —  .  (27) 

dx2 

To  a  reasonable  approximation  the  temperature  of  the  plasma  is  constant, 
T  ~  1  eV;  so  there  is  heat  conduction  from  the  plasma  into  region  B. 
By  analogy,  with  chemical  combustion  waves,  Raizer33  felt  that  heat 
conduction  would  supply  sufficient  energy  to  increase  the  temperature  of 
region  B  to  "combustion".  Careful  calculation  by  Hall,  et.al.18  and 
Thomas,  et.al.28  have  shown  that  this  is  not  the  case. 

(IT)  Jackson  and  Nielsen32  have  investigated  the  heating  of  region 
B  by  radiation  transport  from  the  plasma  of  region  C.  The  opaque  1  eV 
temperature  plasma  emits  strongly  throughout  the  visible  region  of  the 
spectrum  and  also  well  into  the  ultraviolet.  The  dominant  radiation 
heating  contribution  in  region  B  comes  from  this  source  rather  than 
from  the  laser  beam.   Effectively,  the  plasma  serves  as  a  frequency 
converter  for  the  laser  energy,  and  re-emits  it  in  a  spectral  region 
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Fig.  7  Abortion  coefficient  of  a  carbon  dioxide 
laser  in  air  as  a  function  of  temperature. 
From  reference  18. 
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Fie.  8  Thermal  conductivity  of  air  as  a  function  of 
temperature.  From  reference  18. 
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where  the  weaker  plasma  of  region  B  is  already  opaque, 
(U)  The  radiant  flux  term  is  written  as 


CO 

$  =  -c  \   dvK'  (G   CO  -  G  )  (28) 

J      v   vp       V  K      J 


in  the  Jackson  and  Nielsen  model.  Here  G  is  the  irradiance  per  unit 
frequency  from  all  sources  and  G   (T)  is  the  equilibrium  radiative 
energy  flux  density.   In  this  formulation  $  is  a  loss  term  in  region  C 
and  a  source  term  in  region  B.  The  simple  attenuation  approach  of 
Eq.  23  comprises  one  part  of  this  expression.  The  coefficient  K'  is  the 
absorption  coefficient  at  frequency  v  corrected  for  induced  emission, 
see  Zel'dovich  and  Raizer31. 

(U)  This  coupling  of  the  various  portions  of  the  LSC  wave  through 
radiation  transport  makes  the  entire  model  strongly  non-linear.  Still, 
the  model  equations  can  be  solved  analytically  if  three  assumptions  are 
made: 

1)  The  absorption  coefficient  must  be  treated  as  a  step  function 

as  in  Eq.  25.  (Jackson  and  Nielsen  use  0.6  cm  x  with  a  threshold  tempera- 
ture of  1  eV) . 

2)  The  ratio  3  0O/C  (T)  is  set  equal  to  a  constant  f  =  1.0  x  10~3 
gm/cm-sec.  The  enthalpy  of  air  is  shown  in  Figure  9,  together  with  a 
useful  analytic  approximation.  The  relation  C  =  (p)h/3T  )  and  this 
approximation  give  f  w  2  x  10  3  gm/cm-sec. 

3)  The  radiative  transport  term  <S>,  is  written  as  an  exponential 
heating  term  in  front  of  the  plasma  front  and  a  constant  absorption  term 
behind  the  plasma  front. 

(U)  The  system  of  equations  can  then  be  solved  by  a  series  of 
iterations  which  adjust  the  values  of  the  necessary  parameters.  Jackson 
and  Nielsen  find  that  the  most  important  contribution  energy  contributes 
from  radiation  transport  comes  in  the  range  14.5  eV  <^  V  _<  20  eV,  which  is 
well  into  the  uv  region.  Account  must  be  taken  of  the  anisotrophy  of  the 
radiation  field.  The  authors  use  19  frequency  groups  to  establish  the 
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Fig. 9  Enthalpy  of  air  as  a  function  of  temperature. 
From  reference  18. 
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self-consistent  solutions. 

(u)  The  model  preduces  the  correct  trend  of  the  propagation 
velocity,  D,  with  the  laser  irradiance,  3,  but  agreement  in  magnitude 
is  not  satisfactory;  see  Figure  10. 

(U)  Jackson34  has  shown  that  the  presence  of  the  target  surface 
forces  a  streaming  velocity  upon  the  air  in  which  the  LSC  wave  is  propa- 
gating; so  that  the  experimentally  determined  velocity  is  considerably 
larger  than  the  radiation  transport  model  predicts. 

(u)  The  final  transport  term,  £,  the  electron  fluence,  has  appar- 
ently not  been  considered  by  any  investigator.  Clearly  region  C  contains 
a  high  density  of  electrons  and  no  free  electrons  are  present  in  region 
A;  so  there  must  be  a  diffusion  of  electrons  into  region  B  driven  by  a 
concentration  gradient.  The  effect  is  complicated  by  the  fact  that  ions 
must  accompany  the  electrons  to  maintain  charge  neutrality;  so  the  model 
is  one  in  which  a  high  density  plasma  diffuses  into  a  low  density  plasma. 
This  type  of  an  analysis  would  require  an  additional  conservation  law, 
the  number  of  particles,  and  a  kinetic  theory  approach  to  the  spatial 
distribution  of  particles. 

(U)  One  other  important  idea  is  suggested  by  Figure  6,  which 
depicts  both  LSC  and  LSD  waves  produced  by  the  same  computer  code.  With 
the  LSC  wave  irradiation  condition,  the  LSC  wave,  the  plasma  front,  falls 
behind  the  shock  front;  so  the  wave  presents  the  characteristic  double 
structure.   As  the  irradiance  increases,  the  speed  of  the  plasma  front 
increases  until  finally  the  two  front  propagate  together  as  an  LSD  wave. 
Apparently  there  is  a  continuous  variation  from  plasmatron  through  LSC 
behavior  to  LSD  behavior. 

(U)  Jackson35  has  calculated  the  irradiance  at  which  the  LSD  behavior 
should  appear.  He  finds 


12     3  2 
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where  p.  is  the  ambient  air  density  and  hR  is  the  enthalpy  in  region  B. 
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Fig.  10  Plasma  wave  velocity  as  a  function 
of  irradiance. 
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With  reasonable  values  for  the  parameters,  H  =  2  x  107  W/cm2  which  agrees 
quite  well  with  the  experimental  results  of  2.25  x  107  W/cm2  18. 

4.4  LSA  WAVES  IN  TWO-DIMENSIONS 

(U)  The  geometry  of  any  LSA  wave  system  is  clearly  cylindrical  as 
soon  as  it  enters  the  transition  phase.  The  one-dimensional  analyses 
presented  here  are  only  approximations  which  have  meaning  when  radial 
effects  (across  the  wave)  can  be  neglected  relative  to  longitudinal 
effects  (parallel  to  the  direction  of  propagation).  When  ignition,  trans- 
ition, and  propagation  are  under  study  radial  effects  introduce  additional 
loss  mechanisms  which  withdraw  energy  and  electrons  from  the  plasma. 

(U)  Magnetic  confinement  of  the  plasma,  caused  by  its  motion,  will 
be  neglected  in  this  discussion.   It  is  an  additional  consideration  which 
will  tend  to  reduce  the  effect  of  the  radial  loss  mechanisms. 

(u)  The  significance  of  the  radial  direction  is  effectively  deter- 
mined by  the  diameter  of  the  laser  beam  in  the  volume  under  study.  A 
small  diameter  beam  allows  significant  radial  transport  of  electrons, 
heat,  and  radiation;  a  large  diameter  beam  does  not.  Raizer3  has  shown 
that  when  the  amount  of  radial  of  expansion  is  small  the  energy  balance 
equation  should  be  corrected  by  the  replacement  G  ■*  G/(l  +  £/R) ,  where 
R  is  the  radius  of  the  laser  beam,  and  Z   is  the  length  of  the  absorption 
region;  that  is  the  depth  of  the  plasma  in  region  C  of  the  last  section. 
Jackson  and  Nielsen32  have  done  rather  careful  calculations  which  confirm 
that  radial  losses  are  of  secondary  importance  in  the  analysis  of  LSA 
waves  when  the  condition  t  <  <   R  is  satisfied. 
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5.0  SUMMARY 

(U)  This  article  reviews  the  initiation  and  propagation  of  LSC  and 
LSD  waves.   It  concludes  that  prompt  ignition  requires  that  electrons  be 
emitted  from  the  target  surface,  most  probably  at  physical  features 
and/or  inclusions  in  the  surface.  The  electrons  are  followed  by  ions 
extracted  from  the  surface  by  electrostatic  forces.   If  the  irradiation 
is  sufficiently  intense  the  plasma  moves  away  from  the  surface  and  feeds 
upon  the  incident  laser  radiation.  The  plasma  is  accompanied  by  a 
shock  wave.  At  intermediate  irradiance  levels  the  plasma  falls  behind 
the  shock  wave  and  an  LSC  wave  results.   For  irradiance  levels  above 
***  2  x  107  W/cm2  the  plasma  and  shock  wave  move  together  as  an  LSD  wave. 
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